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ABSTRACT 


Observation and energy spectra of electrons and pro- 
tons in islands at the back of the magnetosphere are pre- 
sented. It is shown that typically for electrons J(> 20 kev) 

< 8 x 10 6 /(cm 2 -sec-ster) and for protons J(> 125 kev) < 

10 4 /(cm 2 -sec-ster). The electron spectra for particles 
with 20 kev<E <70 kev can be characterized by e-folding 
energies ranging from 1 5 to 5 kev or by y from 2.6 to 
4.3 for differential power law spectra. For protons with 
energies E > 125 kev the spectra display e -folding energies 
in the range from 20 to 40 kev for y between 7.8 and 4.9 
for differential power law spectra. The electron spectrum 
in the islands shows a characteristic softening with time. 

A comparison is made with electron fluxes observed with 
the Alouette, IMP-1, and Vela satellites. The similarity 
between particles in the islands and those observed in the 


aurora is discussed. 
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ELECTRON AND PROTON FLUXES 


IN THE 

TAIL OF THE MAGNETOSPHERE 
By 

Andrei Konradi 
Goddard Space Flight Center 


INTRODUCTION 

Interesting results concerning high electron fluxes in the tail of the 
magnetosphere were reported by several groups of investigators (Anderson, 
1965; Anderson et al. , 1965; Montgomery et al. , 1965), Measurements 
made directly in the tail by detectors flown on IMP-I and the Vela satel- 
lites indicate that fluxes of electrons with energies above 45 kev are found 
in the tail of the magnetosphere as far out as 31.5 R e . The region in which 
these electrons occur is bounded roughly by the ±20 degree parallel of 
geomagnetic latitude and 130 and 270 degrees ecliptic longitude at 17.7 
R e . The electron fluxes form islands which seem to be time dependent 
phenomena as evidenced by rise times in the flux of a few minutes and 
roughly exponential decays with characteristic times measured in hours. 
They are concentrated mainly close to the trapping region and the fre- 
quency of their occurrence decreases with distance. The spectra of 
electrons observed can be approximated by integral spectra of the form 
E~ n with n about 3.2 to 4.2 between 50 kev and 150 kev. The peak fluxes 
found are a few times 10 6 / (cm 2 - sec- ster ). 
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Observations of electron fluxes outside the trapping region were 
also made in the vicinity of the earth at a height of 1000 km by detectors 
flown on board the Alouette I satellite (McDiarmid and Burrows, 1965). 

l 

High intensity electron fluxes of short duration were observed at high 
latitudes outside the boundary of the outer radiation belt. The intensities 
of these spikes were found at times to reach 10 9 (cm 2 -sec-ster) and the 
electron spectrum observed was softer than the spectrum of electrons 
trapped in the outer belt. The spikes are concentrated in the night side 
of the earth. Typically the latitude width of the spikes is less than 2 
degrees, and some evidence indicates that the spikes are approximately 
aligned with magnetic L- shells. They also occur predominantly during 
times of enhanced magnetic activity. 

The purpose of the present paper is to discuss some additional 
features of the islands: the detection of low energy protons, measure- 

ment of electron and proton spectra and the observation of energy 
dependent time decay of electron intensities. A comparison of these 
features is also made with the work of Anderson et al., Montgomery 
et al. , and McDiarmid et al. A discussion and comparison with electron 
and proton spectra observed in the auroral zone is also given. 

The results reported have been obtained from the ion-electron 
detector flown by Davis and Williamson on Explorer XIY. 
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APPARATUS 

The instrument used in this analysis has been described elsewhere 
(Davis and Williamson 1962, Davis 1965). Here we shall give only a 
brief summary of its operational characteristics. 

The detector is a scintillation counter consisting of a photomulti- 
plier tube on the face of which is deposited a 5 mg/ cm 2 thick layer of 
crystalline ZnS covered by a 1000 A thick layer of Al. 

A wheel driven by a stepping motor introduces varying thicknesses 
of Ni absorber between the collimator and the phosphor, thus permitting 
acquisition of energy information on the incident particles. The wheel 
also carries three Au discs which in some wheel positions scatter the 
incident flux entering through an alternate collimator into the phosphor. 

The detector has a geometric factor of 7.25 x 10' 3 cm 2 - ster for 
direct geometry and 8.13 x 10' 3 cm 2 -ster for scatter geometry. A 
collimator allows particles to enter from a viewing cone with a half 
angle of about 7 degrees. 

There are three modes of operation: 

1. Pulse output from the anode measures protons between 
the energies of 100 kev and 5 Mev. 

2. Normally eighth dynode current measures the total inci- 
dent energy flux from both protons and electrons. 

3. Electrons can, however, be selected preferentially by 
scattering the incident beam of both protons and electrons 
from the Au discs. In this mode of operation the eighth dynode 
current measures the total incident energy flux due to electrons 
only. 
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In all three modes of operation the lower energy cutoff can be 
raised by introducing the varying thicknesses of Ni absorbers mounted 
on the wheel. 

After launch in October, 1962, the detector operated in all three 
modes until the beginning of December, 1962, from which time it failed 
to operate in the first mode, apparently due to the malfunctioning of the 
satellite encoder. Thus no information about the presence of low energy 
protons is available from that time on. 

ANALYSIS 

Due to the directionality of the detector and the spin of the satel- 
lite, it is possible in principle to measure pitch angle distributions of 
particles if the directions of the magnetic field and of the particle 
detector are known as a function of time. While for L < 6 it is possible 
to compute the field direction with good accuracy from a harmonic 
expansion, it becomes necessary to use experimentally obtained results 
on higher L shells and in particular in the tail of the magnetosphere. 

In the case of Explorer XIY the determinations of both the look direc- 
tion of the detector and the direction of the magnetic field were greatly 
complicated by the satellite's precession during the first three months 
of its life. To facilitate data processing, we used only the peak intensi- 
ties recorded during a 5.2 second period, the dwell time for a wheel 
in a given position. For trapped particles within the sunward part of 
the radiation belts where the intensity peaks at 90 degrees to the field 
line this corresponds to looking at pitch angles close to 90 degrees, 
while in other regions this corresponds only to looking in the direction 
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of maximum intensity within the region scanned by the detector during 
the 5.2 second period of observation. The proton spectra presented 
here were obtained by assuming that they can be represented by either 
j(>E) = Ae' E ^ E ° or j(>E) = BE"^ and calculating the E 0 or y from two 
integral intensity measurements with cutoffs at 125 and 168 kev. In 
general the uncertainty in E 0 is ± 4 kev, while the uncertainty inyis 
about ±1. 

The electron spectra were deduced by comparing experimental 
eighth dynode currents with currents computed from 

( 00 

j(E)ES i (E)dE i = 1, 2 10. 

o 

Here G is the geometric factor, j(E) is the assumed differential 
energy spectrum, E is the energy of electrons, and S ; (E) is the sensi- 
tivity curve of the detector for total energy flux with the ith absorber 
in the path of the incident electrons. 

To facilitate comparison with other experimenters, spectra of 
the forms Ae _E ^ E ° and BE'^ were used in the analysis. The para- 
meters E 0 andy were obtained through a weighted least squares fit 
of the calculated currents to the experimental eighth dynode current. 
The experimental uncertainties in the parameters due to the calibration 
and drift in the electronics should be about 10 to 20 percent. Other 
uncertainties in the spectrum are due to short time fluctuations in the 
electron intensity and amount to about 20 percent. 
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The uncertainty in the total electron intensity above some E is 
due to assumptions of different spectral forms and usually runs about 
50 percent. The observations reported here were made in the region 
between 0100 and 0900 local time, 6 and 16 R e , and approximately 
+ 10 and -30 degrees magnetic latitude. 

Two typical electron islands observed by Explorer XIV in the 
night side of the earth at a distance of about 15 R e are shown in Figure 

1. The diagram consists of four curves representing the eighth dynode 
output current for three different absorber thicknesses with cutoffs of 
roughly 20, 30, and 40 kev for electrons. The numbers above the 
graph indicate E 0 (or y given in parentheses) calculated from the 
detector response fon different absorbers. The horizontal bars above 
the curves indicate the range over which the data was averaged to 
obtain the energy spectra. 

There are three features characteristic of electron islands to 
be noted here : 

1. The rise time of the flux is 5 to 10 minutes. 

2. The decay time is of the order of hours. 

3. The decay time is faster for higher energy particles 
implying an overall softening of the electron spectrum. The plot 
is linear in time with the distance (R e ), the geomagnetic latitude 
(GMLAT) and the geomagnetic local time (GMSEP) also indicated. 
Figures 2a and 2b show electron and proton fluxes encountered 

during two outbound passes on 22 and 25 November 1962. The upper 
graph represents the integral intensity of protons with energies above 
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125 kev. The lower graph gives the value of the eighth dynode current 
of the photomultiplier which in this case is a measure of the electron 
intensity for electrons with energies above 20 kev. To simplify the 
diagram only one curve is shown. In Figure 2a the satellite passed 
through the boundary of the trapping region which was located at 
approximately 6 R e in the antisolar direction. 

The close position of the boundary on the night side of the earth 
agrees well with observations made by Anderson et al. (1965) and 
the model presented by Frank (1965). Just outside the boundary 
the satellite passed through a region of low energy electrons with 
E q about 3.7 kev which extended over 1.5 R e . Then after passing 
through several small increases in particle intensity, the first en- 
counter with an island occurred at 0610 UT when the eighth dynode 
current rose by 2.5 orders of magnitude above the background during 
a period of several minutes. The e -folding energy of the electron 
spectrum became 9.1 kev. The intensity of electrons above 20 kev 
corresponding to this spectrum was about 8.4x1 0 6 e /(cm 2 -sec-ster). 
During the encounter with the island the intensity of particles dropped 
and rose two more times. Each of the peaks was lower than the one 
preceding it, indicating an overall decay time of about one hour for 
electrons with E > 20 kev. 

Three more islands were observed by the detector during that 
pass at 0830, 1110 and 1240 UT. It should be noted that in islands 2 
and 3 there is a clear tendency for E q to decrease with time. This 
effect may very well be also present in island 1 but could easily be 
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masked by 10-20 percent temporal fluctuations. To decrease the 
effect of these fluctuations the spectrum in each case was computed 
from an average of several complete spectral readouts. The range 
over which the averaging was done is indicated in the diagram by a 
bar above the curve. Characteristically the e-folding energy observed 
ranged from about 10 to 5 kev. As can be seen from the upper curve, 
protons with energies above 125 kev were observed to coincide with 
the appearance of the electrons. The highest intensities observed 
were about 2 x 1 0 4 p/(cm 2 -sec-ster). The e-folding energy of protons 
ranged from above 20 to 40 kev. The low relative intensity of detected 
protons make it difficult to see whether there is any time dependent 
softening of the spectrum. 

Figure 2b shows another pass when the satellite observed several 
more islands of protons and electrons. Again the same general fea- 
tures observed in Figure 2a are discernible here. Particularly clear 
is the gradual softening of the electron spectrum as a function of time 
during each encounter with an island. The e -folding energies range 
from 16 to 5 kev. Again protons were observed with e-folding energies 
from 30 to 40 kev and intensitie s up to 1 0 4 / (cm 2 - sec-ster). 

The proton fluxes observed do not necessarily coincide exactly 
with the electron islands. Thus the last proton island observed seems 
to end at about 1120 UT while the electron fluxes persist for several 
more hours. This of course may simply be due to the proton flux 
intensity falling below the threshhold of detectability of our detector. 
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Figure 3 shows the appearance of islands on a magnetically quiet 
day, the K p daily sum being 7 0 . On such quiet days and at low geo- 
magnetic latitudes, the boundary of the trapping region sometimes 
extends tol2R e in the antisolar direction, (Serlemitsos, 1965). At least 
one island which started at 1400 UT can be seen in the graph. In 
this instance we can make a direct comparison between an island and 
the edge of the outer belt in the tail region of the magnetosphere which 
the satellite enters at about 1455 UT. It should be noted (1) that the 
total intensity of electrons above ^20 kev is higher in the island than 
in the part of the outer belt adjacent to the island, (2) while initially 
harder, after almost an hour, the spectrum of particles in the island 
softens to the point when it is softer in the island than at the edge of 
the trapping region, (3) while shortly before the satellite passes into 
the trapping region the overall intensity of electrons above approxi- 
mately 20 kev is higher than at the edge of the trapping region, the 
intensitv of electrons above 50 kev is lower. 

The differences in the spectral shape of the electrons found in the 
islands and in the intensities of the electrons indicate that it is unlikely 
that the islands are produced simply by regions containing particles 
becoming detached from the outer belt and drifting into the tail. 

This, of course, does not preclude the possibility that particles 
trapped in the outer belt undergo sudden acceleration and then detach 
from the trapping region. In addition, however, a mechanism must be 
found to explain the energy dependent loss of electrons in order to 
account for a lower intensity of higher energy electrons in part of the 

island than at the edge of the trapping region. 
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Another possible explanation would be that normally closed field 
lines at the edge of the magnetosphere break up into an open configura- 
tion, allowing the trapped particles to be lost and bringing in particles 
accelerated at a location as yet unknown. 

A test of these hypotheses will have to await a more detailed 
analysis involving correlations with data from the onboard magnetom- 
eter. 

DISCUSSION AND CONCLUSIONS 

The following are some of the characteristic features of the 
particle fluxes in the tail of the magnetosphere: 

1. Electron fluxes have rise times of the order of 5 minutes. 
The decay times, however, range from 30 minutes to several 
hours depending on the energy of electrons measured, with the 
higher energy electrons decaying faster than the lower. 

2. The E 0 of electrons is of the order of 1 0 kev (y of the 
order of 4). 

3. The electron spectrum softens with time with the e- 
folding energy changing from about 1 1 to 6 kev over a period of 
one or two hours suggesting an energy dependent decay mechanism. 

4. Associated with the electron fluxes are fluxes of protons. 
The E q for protons range from 20 to 40 kev. There is no clearcut 
change with time in the proton spectrum observed. Also, at 
times while there exists a clearly discernible electron flux, the 
proton flux apparently falls below the threshhold of detectability 
which is 10 ? p/(cm 2 -sec-ster). 
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5. Calculations show that maximum intensities of electrons 


with E > 20 kev may reach about 8xl0 6 e/(cm 2 -sec-ster). 

The electron fluxes observed in the tail of the magnetosphere 
have a strong resemblance to the islands discussed by Anderson 
(1965). In particular the rise times of several minutes and decay 
times of hours seem to be identical. The intensities quoted by 
Anderson, j(E)< a few times 10 7 /(cm 2 -sec), are quite comparable to 
those reported here since if we raise the cutoff energy in our case 
from 20 to 45 kev so as to compare with the Geiger counter used by 
Anderson, and calculate the omnidirectional flux we get~10 7 e/(cm 2 - 
sec) which compares well with the upper limit of intensity observed 
by Anderson. We shall henceforth assume that the electron islands 
observed by Anderson and those described here are identical. 

Comparison of our fluxes with those observed by Montgomery et 
al. (1965) on the Vela satellites also indicate similarities. The 

fluxes observed by Montgomery et al. are confined to approximately 
±20 degrees magnetic latitude. A preliminary survey of our data 
also indicates that we seldom see electrons at magnetic latitude 
greater than -25 degrees. (The maximum excursion of the orbit in the 
northern hemisphere is about .-^degrees geomagnetic latitude for 
distances >10 R e . Thus we cannot sample the northern latitude limit). 
The total electron intensities reported range typically from 10 4 to 
10 5 e/(cm 2 -sec). Again since the cutoff energy of the detector used 
was 50 kev, the values are quite comparable to those observed by us. 
Comparison of the electron energy spectra shows that typical values 
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of n obtained by fitting an integral energy spectrum of the form 
j(>E)~E' n to the observed electron spectra gives values of n = 3.2, 

3.4 and 4.2. Presumably this fit was done in the energy range between 
50 and 150 kev. This would correspond to y (for a differential energy 
spectrum where n = y +1) of 4.2, 4.4, and 5.2. On the other hand if 
only a two point spectrum is used at 50 and 70 kev, which is closer 
to the present measurements, one obtains values of y (for a differential 
energy spectrum) equal to 3.24, 3.8, and 4.8 which are quite close to 
those reported here. While not stated explicitly, a look at Figure 2 
of Montgomery et al. (1965) shows a gradual softening of the 
spectrum with time comparable to the softening reported here. This 
can also be seen from the change in the exponent n or y used to repre- 
sent the spectrum. The division into early, main, and late phases used 
by Montgomery et al. can also be applied to some of the electron fluxes 
observed here, e.g. Figure 2 between 0640 and 0930 UT. No compari- 
son between time fluctuation can be made since the resolution time 
under data processing conditions reported here is 83 seconds. 

McDiarmid and Burrows (1965) reported observations of elec- 
tron spikes made in the night side at a height of 1000 km beyond 
the outer radiation belt approximately between the invariant latitudes 
A = 71 and 82. A two point fit of electron intensities above 40 kev 
and 250 kev produces an exponent in a differential power law spec- 
trum ranging from > 3.6 to> 7. 

The measurements reported here show a range of y between 
approximately 2.6 and 4.6 which are somewhat lower than those of 
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McDiarmid and Burrows. On the other hand in our case the spectral 
fit was made between 10 and 50 kev. The reported intensities range 
from ~ 2. 1 0 4 e / (cm 2 -sec-ster) at A = 82 degrees to ~ 8x1 0 8 e / (cm 2 - 
sec-ster) at A = 70 degrees. Considering that the cutoff energy of 
the detector is 40 kev, fluxes of 10 s /(cm 2 -sec-ster) are more than an 
order of magnitude higher than those observed here. Since the detec- 
tor used in the present study has an aperture of ± 1 1 degrees it is 
possible that it would average out a highly anisotropic flux streaming 
down the field line. (It should be remembered a mirroring particle 
at 1000 km at A =70 corresponds to a particle with a pitch angle of 
only 1.4 degrees at the equator for the earth's dipole field). The 
effect of this would be to decrease the magnitude of the apparent 
observed flux near the equator. There is some evidence reported 
of a possible time decay in the intensity of observed particles based 
on the difference in intensities observed at the same L on the ascend- 
ing and descending leg of the orbit. However, since the decay is 
coupled with a change in local time, a meaningful comparison is 
hard to make. We feel therefore the need of more observations before 
we can attempt to identify the spikes observed by McDiarmid and 
Burrows with the electron intensities observed here. 

Another interesting comparison can be made with the intensities 
and spectra of electrons and protons precipitated in the auroral zone. 

An excellent summary of auroral measurements done up to 1964 
was made by Hultqvist (1964). The auroral observations can be 
divided between those done from balloons, rockets and polar orbit- 
ing satellites. 
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The range of energies investigated covers a region extending 
from 3 to 250 kev for electrons and 80 to 800 kev for protons. The 
spectra observed are highly variable in time and in general cannot be 
approximated by a simple exponential or power law over the whole 
range. Thus it seems reasonable for purposes of comparison to 
restrict our attention only to measurements over energy ranges which 
substantially overlap the range reported here and where either a 
particle energy spectrum is reported or can be obtained from other 
data cited. 

Table 1, partially adapted from Hultqvist, presents a summary of 
such spectral measurements. The outstanding feature of the auroral 
measurements is the variation in the shape of the spectra and the 
intensity of precipitated fluxes. The spectral form changes both with 
energy and with time. Virtually no other detector measured particles 
over the same range of energies as our detector. One can therefore 
not compare spectral parameters {y and E q ) between two measure- 
ments and expect full agreement or disagreement. For example, a 
spectrum that is close to exponential with E q ~ 1 0 kev and which 
was determined from an integral measurement at say 10, 40, and 
90 kev might be approximated with two different power law spectra 
in the range 10 to 40 kev and 40 to 90 kev - namely y =1.17 and y - 
5.17 which differ considerably. Under these circumstances a com- 
parison can be made only subject to the following limitations: 

1. Only measurements done over ranges of energy sub- 
stantially overlapping our range can be considered (Table 1). 
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2. Measurements producing either E 0 or y similar to those 
observed by us shall be considered in agreement in the hope that the 
observed spectrum favors either an exponential or a power law. Sub- 
ject to these conditions good agreement is observed with auroral measure- 
ments done by Mcllwain (I960), Bhavsar (1962), Anderson and DeWitt 
(1962), Krasovskii et al. (1962), O'Brien and Laughlin (1962), O'Brien 
et al. (1962), Stilwell (1963), Sharp et al. (1964), McDiarmid and Budzinski 
(1964). In the case of Davis et al. (I960) the spectra both for protons and 
electrons are somewhat harder than those observed here. If, however, 
an exponential form is assumed for Davis' electron spectrum, E 0 comes 
out to be approximately 10 kev which is in the range of spectra observed 
by us (private communications). Harder spectra were also observed by 
McDiarmid et al. (1961), Anderson and Enemark (I960), Sharpet al. 

(1964) and Evans (1965). 

Hultqvist (1964) has used data obtained by O'Brien (1964), O'Brien 
and Taylor (1964) and McDiarmid et al. (1963) to calculate e-folding 
energies from average measurements obtained from Injun III and Alouette 
over the auroral zone. The results are E 0 ~ 5.7 kev (1 kev < E < 40 
kev) and E Q 41 or 31 kev (40 kev < E < 250 kev). As pointed out by 
Hultqvist one should keep in mind the questionable meaning of averages 
and the rather crude approximations made in the calculations. Nonethe- 
less if one takes the numbers at their face value, the results obtained 
here are in good agreement with them. This survey of auroral data 
shows that there exists a great deal of similarity between the electrons 
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and protons observed in the aurora and in the islands. It is tempting, 
therefore, to identify them with each other. 

Observations made by Ness (1965) and Anderson (1965) indicate 
that the islands occur on apparently open field lines; however, O'Brien's 
(1964) observation of 1.5 Mev electrons trapped during auroral precipi- 
tation of lower energy electrons leads to the belief that aurora occurs 
on closed field lines. Therefore, it would be necessary to invoke a break- 
down of the second adiabatic invariant to permit particles to cross mag- 
netic shells. 

However, we feel that lacking more information about pitch 
angle distributions, spatial distributions, time variations, and the 
detailed shape of the magnetic field in the tail of the magnetosphere 
one cannot make any statement at this time about the possible casual 
relationship between fluxes in the magnetospheric tail and the aurora. 

We would like to suggest, however, that the similarity of the spectra 
and intensities of both electrons and protons in the aurora and the 
islands makes it probable that a common accelerating mechanism 
is responsible for their production. 
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Summary of Spectral Measurements 
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LIST OF FIGURE CAPTIONS 


Figure 1— Appearance of two typical electron islands. The four curves 
correspond to the ion-electron detector response with four 
different absorbers in the path of the incident beam. The 
numbers above the curves give the value of E 0 (or y in 
parenthesis). 

Figure 2 — Outbound passes of (a) 21 November and (b) 25 November 
1962. In both cases, the upper curve gives the maximum 
intensity of protons with E > 125 kev, while the lower curve 
corresponds to the maximum ion-electron detector response 
to electrons with ¥ 20 kev. The numbers above the curves 
indicate the value of E Q (or y in parentheses). 

Figure 3 — Appearance of an electron island shortly before the satellite 
entered the trapping region (outer belt). The four curves 
correspond to the ion-electron detector response with four 
different absorbers in the path of the incident beam. The 
numbers above the curves give the value of E 0 (or y in 
parentheses). 
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